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1 Introduction 

The production of 77 mesons and their collisions with nuclei have been studied experimen- 
tally and theoretically with increasing interest during the last years. To a large extent 
this is motivated by the fundamental problem of charge-symmetry breaking of strong 
interaction pQ. Another relevant questions concerns the nature of Sii(1535)-resonance j2] 
and the possible formation of ^-nucleus quasi-bound states It is worth mentioning 
that according to [3] the mean free path of 77 mesons in a nuclear medium is about 2 fm, 
i.e., less than the size of a typical nucleus. A necessary condition for the existence of 
77-nuclei, hence, appears to be satisfied. 

For the calculation of these states various model treatments were employed, among 
them the optical potential method the Green's function method |7], the modified 

multiple scattering theory jHj and few-body calculations [H] - |Hj. The predictions con- 
cerning the possibility of 77-mesic nucleus formation are very diverse. One obvious reason 
for such a diversity is the poor knowledge of the i]N forces. Another reason comes from 
the differences among the employed approximations some of which might be faulty in 
view of the resonant character of the 7]N dynamics and the delicacy of the quasi-bound 
state problem. 

In the present paper we treat the 77-deuteron system on the basis of the exact few- 
body Alt-Grassberger-Sandhas (AGS) equations. The Faddeev-type coupling of these 
equations guarantees uniqueness of their solutions. Moreover, as equations for the elastic 
and rearrangement operators they are well-defined in momentum space, providing thus 
the desired scattering amplitudes in a most direct and technically reliable manner. The 
advantage of working with coupled equations involving the elastic and rearrangement 



operators of the final state, is not only suggested by questions of uniqueness, but also by 
the relevance of rescattering effects which were found to give a significant contribution to 
the corresponding amplitude [12]. 

2 AGS formalism 

In terms of the AGS transition operator U\\ the 7]d elastic scattering amplitude is repre- 
sented as 

/(pi.pi;*) = -(27r) 2 M 1 (p / 1 ;V ci |[/ii(^)ipi;^) (1) 

with the on-energy-shell conditions \p[\ = |pi| and z = p\/2M\ + Ed with Ed being the 
deuteron energy. Here the subscript 1 labels the r](NN) partition and the ^-deuteron 
channel. The transition operator U\\ obeys the system of AGS equations 

3 

Upa(z) = (1 - 5 Pa )G \z) + ]T(1 - 6^)T^z)G (z)U ia (z), (2) 

7=1 

with Gq(z) being the free resolvent (Green's operator) of the three particles involved. 
This set of equations couples all 3 x 3 elastic and rearrangement operators U aa and 
Up a . Here each of the subscripts runs through the values 1, 2 and 3, indicating the 
two-fragment partitions (1,23), (2,31) and (3,12) respectively. Therefore, Un describes 
the elastic transition 1(23) — > 1(23), while U21 represents the rearrangement process 
1(23) — > 2(13). These subscripts are also used in the complementary notation to label 
the two-body T-operator t a (z) of the (^7) pair, for instance t\{z) — ijvjv(z). It should 
be noticed, however, that it is not this genuine two-body operator which enters the AGS 
equations, but the operator 

T a {z) =t a {z-ql/2M a ), (3) 

which is to be understood as the two-body operator embedded in the three-body space, 
with the relative kinetic energy operator q£ l /2M a of particle a being subtracted from the 
total energy variable z. Considered in momentum space, (j3J) thus reads 

(p'a^a\ T a(z)\p a ,q a } = 5(q! a - q a ) (p' a \t a (z - q 2 j2M a ) |p a ) , (4) 

where p a and q Q are the Jacobi momenta of the pair a and the spectator a respectively, 
and M a is the corresponding reduced mass. 

For both T vN and T NN we used one-term separable forms 

T a {z) = \Xa)r a (z)(xa\ ■ (5) 

For the NN subsystem Eq. (jSJ implies that the asymptotic wave function is related to 
the form-factor \xi) according to 

\pi,^ d ) = G (z)\xi)\pi) ■ (6) 
Due to (0) and (0) the scattering amplitude (JTJ) can be rewritten as 

/(p' a ,Pi;z) = -(27r) 2 M 1 (p' 1 |X 11 W|pi) , (7) 
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where the operators Xp a , defined as 

Xf3a(z) = (xp\G (z)Up a (z)G (z)\x ol ) 
obey the system of equations 



3 ( Pj \ 



with 

Zf3a(z) = (1 - Sp a ){x/3\G (z)\Xa) ■ 

The identity of the nucleons implies that X31 = X21, T3 = t 2 , and Z31 = Z 2 \, which 
reduces the system (JSJ) to two coupled equations 



X n (z) = 2Z 12 (z)r 2 (z - X 21 (z) 



X 21 (z) = Z a W + ^21(^1 (2 - J-] Xu( 2 ) + Z 23 (.)r 2 f 2 - ^ I 21 (z) 



(9) 



Eventually, after making the S'-wave projection of the matrix elements (p^X^lp^) and 
(p'JZpalpa), we end up with one-dimensional integral equations which can be solved 
numerically. 



3 Two— body T— matrices 

The S'-wave nucleon-nucleon separable potential is adopted from Ref. [T3] with its param- 
eters slightly modified to be consistent with more recent NN data (see Ref. |1 lj ) . The 
7/-nucleon T-matrix is taken in the form 

t vN (p\p-z) = (p' 2 + a 2 )- 1 {z _ Eq \ - r/2) (P 2 + a 2 )- 1 (10) 



consisting of two vertex functions and the Sn-propagator in between [10]. It corresponds 
to the process rjN — ► Su — > r)N which at low energies is dominant. The range parameter 
a = 3.316 fm -1 was determined in Ref. while Eq and T are the parameters of the 
Sn resonance [To] . 

E = 1535 MeV - (m N + m v ) , T = 150 MeV . 
The strength parameter A is chosen to reproduce the ^-nucleon scattering length a vN , 

a\E - l T/2) 

the imaginary part of which accounts for the flux losses into the nN channel. 
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The two-body scattering length a n N is not accurately known. Different analyses JH] 
provided for the values in the range 

0.27 fm < Re a vN < 0.98 fm , 0.19 fm < lma vN < 0.37 fm . (12) 

Recently, however, most of the authors agreed that Im a v ^ is around 0.3 fm. But for Re 
a vN the estimates are still very different (compare, for example, Refs. [T7] and [IE])- Most 
of the results presented in this paper are, therefore, obtained using Im a V N = 0.3 fm and 
several values of Re a V N within the above interval. 

We solved Eqs. © for r\d collision energies varying from zero (r/ci-threshold, z = E^) 
up to 22 MeV by replacing the integrals by Gaussian sums. As is well known (see, for 
example, JH|), the kernels of these equations, when expressed in momentum represen- 
tation, have moving logarithmic singularities for z > 0. In the numerical procedure, we 
handle it with the method suggested in Ref.|2Q|. The main idea of this method consists in 
interpolating the unknown solutions (in the area covering the singular points) by certain 
polynomials and subsequent analytic integration of the singular part of the kernels. 



4 Scattering results 

The results of our calculations are presented in Figs. EE) and in Table ^ In Fig. ^ the 
energy dependence of the r\d phase-shifts for five different choices of Re is shown. 
The curves correspond (starting from the lowest one) to Re a^N = 0.55 fm, 0.65 fm, 0.725 
fm, 0.75 fm, and 0.85 fm. The larger this value, the stronger is the 77N attraction. The 
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Figure 1: Real part of the 7/-deuteron phase- 
shift function of the collision energy. 



Figure 2: The critical values of a V N (filled 
circles) to the right of which the rjd- 
system can be bound. 



change in the character of these curves, hence, reflects the growth of the attractive force 
between the t] meson and the nucleon. The lower three curves for Re S v d corresponding to 
the smaller values of Re a V N start from zero, the two curves corresponding to the strong 
attraction start from tt. According to Levinson's theorem, the phase shift at threshold 
energy is equal to the number of bound states n times ir. We found that the transition 
from the lower family of the curves to the upper one happens at the critical value Re 



4 




Figure 3: Total cross-section for elastic t]d Figure 4: Argand plot for the rjd elastic 
scattering as a function of collision energy. scattering amplitude. 



a>riN = 0.733 fm. Therefore, the r/N force, which generates Im a V N = 0.3 fm and Re 
a V N > 0.733 fm, is sufficiently attractive to bind r\ inside the deuteron. 

For three other choices of Im a n N within the uncertainty interval, namely, 0.20 fm, 0.25 
fm, and 0.35 fm, the corresponding critical values of Re a vN turned out to be 0.788 fm, 
0.761 fm, and 0.698 fm. In the complex a^-plane (see Fig. |2J) the corresponding points 
form a curve separating the uncertainty area (given by formulae (j!2j) . dashed rectangle) 
into two parts. If a V N is to the right of this curve, the strength of 77N attraction is sufficient 
for rjd bound state formation. 



Re (f m ) 


E% (MeV) 


IV (MeV) 


0.55 


8.24 


9.15 


0.65 


7.46 


8.45 


0.675 


7.14 


7.61 


0.70 


6.79 


6.90 


0.725 


6.41 


6.31 


0.75 


6.01 


5.87 


0.85 


4.39 


5.79 


0.90 


3.73 


6.81 



Table 1: Energy and width of the 7]d resonance for various choices of Re a V N. 



In Fig. 01 we present the result of our calculations of the total cross-section (integrated 
over the angles) for elastic 77-deuteron scattering as a function of collision energy. The 
five curves correspond (starting from the lowest one) to Re = 0.55 fm, 0.65 fm, 0.725 
fm, 0.75 fm, and 0.85 fm. The dashed line indicates the deuteron break-up threshold. 
The peaks in the energy dependence of the total elastic cross-section indicate that a 
resonance appears in the rjd-system. Of course, not every maximum of the cross-section 
is a resonance, but we plotted the Argand plots also (Fig. @| for the rjd elastic scattering 
amplitude in the energy interval from to 22 MeV. The five curves correspond (from right 
to left) to Re a V N = 0.55 fm, 0.65 fm, 0.725 fm, 0.75 fm, and 0.85 fm. When the energy 
increases the corresponding points move anticlockwise. So that the Argand plots prove 
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that the maxima we found are resonances. Their positions and widths for various choices 
of Re a V N are given in Table [TJ It should be noted that, while the resonance energy is 
determined in our calculations exactly (as the maximum of the function sin 2 Re S v d), the 
corresponding width is obtained by fitting the cross-section with a Breit-Wigner curve. 
Therefore, the values of T v d given in Table^should be considered only as rough estimates. 

The resonant behavior of rjd elastic scattering should be seen in various processes 
involving ^d-system in their final states, such as 'yd — > rjd and np — > r]d. Indeed, the 
corresponding amplitudes (VwlClVv) involve the rjd wave function ip out which, in the 
vicinity of the resonance, strongly depends on the total energy. Its resonant growth 
at short distances may enhance the transition probability. To check this suggestion we 
performed the calculations of the jd — > rjd reaction. 



5 Photoproduction of ?7-mesons 

Theoretical analysis of (7, /^-reactions on nuclei is hampered by the three major prob- 
lems: the unknown off-shell behavior of the two-body 7iV — > rjN amplitude, inaccuracies 
in the description of the nuclear target as a many-body system, and rescattering effects 
in the final state. The simplest is, of course, the process of coherent rj photoproduc- 
tion on deuteron. There are many theoretical studies devoted to (7, rj) reactions on 
deuteron. Early attempts to go beyond a simple impulse approximation led to very differ- 
ent conclusions [21] - J2H] as do more recent approaches based on the effective two-body 
formulations [21], [2H]- Moreover, the experimental cross-section j2H] of the reaction 

7 + d -> rj + d (13) 

in the near-threshold region is far above these theoretical predictions. Therefore, a reliable 
description of rj photoproduction on deuteron on the basis of exact equations is desirable. 



6 Formalism 

To consider reaction ([13)1 . we employ the exact Alt-Grassberger-Sandhas formalism mod- 
ified to include the electromagnetic interaction. A photon can be introduced into this 
formalism by considering the rjN and 'yN states as two different channels of the same sys- 
tem. This means that we should replace the T-operator t V N by the 2x2 matrix. It is clear, 
that such replacements of the kernels of integral equations ([2]) lead to the corresponding 
solutions having a similar matrix form 



t 77 t 7 M „ / W2J w2 v 

/■rn YXTf. 
a/3 vv ap 



Here t 77 describes the Compton scattering, the photoproduction process, and t m the 
elastic rjN scattering 

1 Another method of treatment electromagnetic process of the type 113(1 in the frame of AGS equations 
was used in Ref.[27|. 
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It is technically more convenient to consider the reaction of 7/-photoabsorption, which 
is inverse to reaction (JT3J). Then the photoproduction cross-section can be obtained by ap- 
plying the detailed balance principle. We are interested in the coherent process, therefore 
we need the amplitude obeying the equation 

( W$ Wfl \ ( T 2 77 \ _ / W£> W? v 



(15) 



In the first order on electromagnetic interaction: 

Wg « T^G (z)W^ + TrG (z)W2? (16) 

and we can see, that in above approximation the transition operator can be found 
not from an equation, but from the expression (jlfij) . It is physically clear since in the 
first order on electromagnetic interaction the only contributions are from the //-meson 
rescattering which is described by the W^x and W^i operators. 

It was experimentally proven [2E| that at low energies the reaction 'jN — > rjN mainly 
goes via formation of the Sn-resonance, which means that t iri in the near-threshold region 
can be written in a separable form similar to (JHJ). To construct such separable T-matrix, 
we used the results of Ref. |2S] where t lv was obtained as an element of a multi-channel 
T-matrix which simultaneously describes data for the processes 

tt + N it + N, vr + N -> r] + N, 
7 + iV -> vr + N, j + N rj + N 

on the energy shell in the ^n-channel. For our calculations, we extended this T-matrix 
off the energy shell, 

QtV>n E ) = n:{E) (17) 

jj K z _|_ p> Q; + p 

using the Yamaguchi form-factors which become unit on the energy shell. Here k stands 
for some unknown parameter. It is known, that t iri is different for neutron and proton, 
we assumed that they have the same functional form and differ by a constant factor 
ty? = AtJ . A multipole analysis [SU] gives for this factor the following estimate: A = 
-0.84 ±0.15. 



7 Photoproduction results 

As was expected, our calculations revealed very strong final state interaction in the re- 
action (|13|). A comparison of the corresponding cross-sections obtained by solving the 
AGS equations and by using the Impulse Approximation (IA) is given in Fig. where the 
IA-results are multiplied by 10. Besides the fact that the IA-curve is generally an order of 
magnitude lower, it does not show a resonant enhancement which is clearly seen when all 
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Figure 5: Total cross-section, calculated 
within a rigorous few-body theory (AGS) 
and Impulse Approximation (IA). 



Figure 6: Differential cross- section 
(9^ m = 90°), calculated with different 
choices of the ratio A. 



the rescattering and re-arrangement processes are taken into account. In this connection, 
it should be noted that experimental data for the reaction (fT3J). given in Figs. El El 
show a pronounced enhancement of the differential cross-section at low energies. 

In order to examine a dependence of our calculations on the choice of the parameters 
of the T-matrices t vv and f 171 , we did variations of A = t'J^/tJ n , Re cl v n and k within 
the corresponding uncertainty intervals. One of the most important parameters of the 
theory is the ratio of the photoproduction amplitudes for neutron and proton (A). Six 
curves corresponding to different choices of A are depicted in Fig. El These curves were 
calculated with a n N = (0.75,0.30) fm and k = a = 3.316 fin -1 . The experimental data 
are taken from Ref. j26j . 

In Fig. [?1 we present the result of our calculations for five different choices of Re 
a V N, namely, 0.55 fm, 0.65 fm, 0.725 fm, 0.75 fm, and 0.85 fm. This sequence of Re 
corresponds to the upward sequence of the curves in the near-threshold region. The 
parameters are: A = —0.75 and k = a = 3.316 fm -1 . 



do/d£2 (nb/sr) 
20 



do/d£2 (nb/sr) 
20 




Figure 7: Differential cross-section (0^ m = Figure 8: Differential cross-section with 
90°) with different choices of Re a V N- different choices of parameter k. 

In Fig. |H1 the results of our calculations with three choices of the range parameter k 
(2 fm -1 , 3 fm -1 , and 5 fm -1 ) of the electromagnetic vertex 7^ — > Su are given. These 
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curves correspond to a n N = (0.75,0.30) fm and A = —0.75. A comparison of the curves 
depicted in Figs. |H1 El and |H1 with the corresponding experimental data shows that no 
agreement with the data can be reached unless the ratio A is greater than —0.80. 
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Figure 9: Shift of the resonant peak of the 
total cross-section due to the rjN repulsion. 



Under all variations of the parameters, however, the resonant peak remains about 15 
MeV to the left of the experimental peak. Since this peak is due to the resonant final state 
interaction between the T] meson and deuteron, we may expect that it can be shifted to the 
right by introducing a repulsion into the rjN interaction. To introduce an rjN repulsion 
which preserves the separable form of the corresponding T-matrix, we used the method 
suggested in Ref. ^3] where a separable nucleon-nucleon T-matrix includes an energy 
dependent factor, b(E) = — tanh(l — E/E c ). This factor causes the NN phase-shift to 
change sign at the energy E c = 0.816 fm^ 1 , which is equivalent to presence of an NN 
repulsion. Since the purpose of our numerical experiment was to check if an rjN repulsion 
could shift the peak to the right and there is no information about such repulsion, we used 
the same function b(E) and did variations of E c , namely E c /3, E c /2, E c , 2E C , and 3E C . 
The corresponding curves are shown in Fig. El where the larger E c the lower is the curve. 
All the curves depicted in Fig. |U] were calculated with a V N = (0.75,0.30) fm, A = —0.85 
and k — ex. — 3.316 fm -1 . 

Therefore, comparison of our calculations with the experimental data suggests that 
A > —0.80, Re &r]N > 0.75 fm, and the rjN interaction is likely to be repulsive at 
short distances. Our calculations are highly sensitive to the ratio of the photoproduction 
amplitudes for neutron and proton (A) and depends rather weakly on the parameter k. 
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